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We describe the electro-thermal transport in metallic carbon nanotubes (m-CNTs) by a semi- 
classical approach that takes into account the high-field dynamical interdependence between charge 
carrier and phonon populations. Our model is based on the self- consistent solution of the Boltzmann 
transport equation and the heat equation mediated by a phonon rate equation that accounts for 
the onset of non-equilibrium (optical) phonons in the high-field regime. Given the metallic nature 
of the nanostructures, a key ingredient of the model is the assumption of local thermalization of 
charge carriers. Our theory remarkably reproduces the room temperature electrical characteristics 
of m-CNTs on substrate and free standing (suspended), shedding light on charge- heat transport in 
these one dimensional nanostructures. In particular, the negative differential resistance observed 
in suspended m-CNTs under electric stress is attributed to inhomogeneous field profile induced by 
self-heating rather than the presence of hot phonons. 

PACS numbers: 73.63.Fg, 73.23.-b, 65.80.+n 



I. INTRODUCTION 

Since their discovery 1 , carbon nanotubes (CNTs) have 
been subject of intense research for both scientific inter- 
est and technological purpose. Their thermal, mechanical 
and electrical properties 2 make them prominent candi- 
dates for broad range of applications in nanotechnology 3 . 
The electronic structure of these one-dimensional (ID) 
macromolecules exhibit semiconducting or metallic (m) 
features depending on their chirality 4 . The latter have 
been proposed for interconnect applications because they 
can support currents densities larger than copper with- 
out suffering electro-migration 5,6 . Recent experiments 
on individual m-CNTs in the high field regime^ ^ 10 ! 11 
have revealed strong nonlinear transport characteristics, 
such as current saturation and negative differential re- 
sistance, which have been attributed to the emergence 
of optical phonon (OP) populations out of equilibrium 
(hot phonons) 1Q i 12 . Additionally, the large electrical 
power dissipated in the high bias regime have fostered 
further studies of the electro-thermal transport in these 
structure o 13 i 14 i 15 i 16 i 17 i 18 . As a consequence numbers of 
theoretical and experimental works have attempted to 
quantify both the electron-phonon (e-p) and phonon- 
phonon (p-p) coupling interaction in m-CNT a 19 i 2Q i 21 i 22 

Prominent among the theoretical approaches for in- 
vestigating high-field transport in m-CNTs at room 
temperature are the Boltzmann transport equation 
(BTE) 2 ^, Monte Carlo simulations^, Landauer- 
Buttiker formalism 10 based on the Fermi golden's rulei 2 - 
that successfully described transport in nanoscale semi- 
conductors. However one of the major drawbacks of these 
approaches is the absence of electron-electron (e-e) inter- 
action, despite the metallic character of the system, and 
the contribution of the charge carriers to thermal trans- 
port. In this work, which is applicable to any metallic ID 



system, we describe the system dynamics as the interplay 
between carrier, acoustic phonon (AP) and OP popula- 
tions. Our self-consistent approach describing electro- 
thermal transport accounts for the onset of phonon and 
carrier populations out of equilibrium and provides a 
comprehensive picture of the carrier flow and heat ex- 
change in the nanostructures. Our model is in excellent 
agreement with experiments on individual m-CNTs in 
both substrate supported (SS) and free standing (FS) 
configuration at room temperature. 



II. MODEL 

We describe the high- field transport in m-CNTs by ac- 
counting for both electric and thermal transport through 
the ID structure. While the former originates from 
the charge carrier populations only, both electron and 
phonon populations contribute to heat transfer in a con- 
jugated manner. 



A. Electric and energy flow 

In small diameter m-CNTs, the higher conduction and 
lower valence subbands lie far above or below the Fermi 
level so that the charge flow arises solely from the contri- 
bution of the conducting bands crossing the Fermi level 
(here taken as the zero-energy level), which are well de- 
scribed by a linear E — k dispersion: 



E±(k) = ±Hv F k 



(1) 



where vp is the Fermi velocity (vp ~ 8 x 10 T cm/sec in 
CNTs) 4 and the + (— ) sign corresponds to the forward 
(backward) carrier branches. Each of these branches has 
a degeneracy 2g c (where the factor of 2 indicates the spin 
degeneracy and g c = 2 for m-CNTs). Due to the efficient 
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FIG. 1: Left: Schematics of the forward and backward quasi- 
Fermi level profiles as a function of the distance along the 
m-CNT. Center: Forward and backward distribution func- 
tions at position x = xq. Right: Band structure at x = xq. 
Contributions of upper and lower subbands are neglected as 
they lie far from the active scattering area (shaded region). 



intrabranch e-e interaction 25 carrier populations are ther- 
malized and described by Fermi-Dirac distribution func- 
tions. In the presence of an external electric field F, the 
local quasi-Fermi levels of forward and backward popu- 
lations separate ^ and the carrier distributions 
in each branch are not necessarily in local thermal equi- 
librium with each other (T+(x) 7^ T_(x)), as illustrated 
in Fig. [TJ The electric current in m-CNTs reads 23 : 



1 g c G 



(2) 



where e is the electron charge and Go = e 2 /(7rh) is the 
quantum conductance. Here we point out that the quasi- 
Fermi level difference does not necessarily coincide with 
the drain-source voltage, i.e. |/i + — < |eV^ s |, where 
the equality only holds in the absence of back-scattering. 
In other words carriers injected from the contacts are 
immediately thermalized by the effective e-e interaction. 
Hence, even if eVd s is larger than the separation between 
subbands, contribution of conducting subbands can be 
neglected provided that the local quasi-Fermi levels in 
each branch lie far away from the upper or lower sub- 
band (\E S — fJL v \ ^> ksTrj), as illustrated in Fig. [TJ We 
also recognize that quasi-Fermi levels fi± (x) are position- 
dependent, but their difference remains constant for cur- 
rent conservation^ 3 -. 

In this metallic system, the distribution function 
f±(E,x) for each branch obeys the stationary BTE ex- 
pressed as: 



± v F [d x f±(E, x) + eF8 E f ± (E, x)} = d t f± 



coll 



(3) 



after using the linear dispersion (Eq. [T]). Here, x de- 
notes the position along the m-CNT and the right hand 
side of the equation is the collision integral. The local 
temperatures and quasi-Fermi levels of the thermalized 
distributions can be determined using the method of mo- 
ments (from here on we omit the variable x for sake of 
brevity) . We have recently shown that the electric field F 



and carrier temperature T± profiles fulfill the relations 25 : 



eF 



, ^ h T coll 



and 



g c nk%T ± IF M 



(4) 



(5) 



respectively. The absence of thermal gradients in Eq. [H 
asserts the lack of thermoelectric power in systems with 
linear dispersion 25 . In Eqs.[4]and[5l ks is the Boltzmann 
constant and 



n^coll 
J~ m,± 



9c 



nhvj 



f 

j — c 



d t f± 



coll 



(E - fi) m dE (6) 



is the m-moment of the collision integral accounting for 
carrier scattering. We express the boundary conditions 
for T± (x) in terms of the transmission coefficient through 
the contacts t, such that: 



Tl(TL/2) = tTl + (1 - t)Tl{TL/2). 



(7) 



where To is the temperature at the contacts. The value 
of t ranges between and 1. The latter corresponds to 
perfect contacts that inject carriers from the leads with- 
out any reflection. After obtaining the F-profile (Eq. 
the net drain source bias is computed as: 



V ds = IR C 



,L/2 

/ F(x)dx 

J-L/2 



(8) 



where R c is the contact resistance which in the case of 
perfect contacts is e 2 /(7rh) 26 . 



B. Electron-phonon interaction 

The main contribution to the collision integrals (Eq. [6j) 
in m-CNTs, other than e-e interaction that thermalizes 
the branch distributions, is scattering with lattice vibra- 
tions (phonons). The collision integrals read: 



dtfn 



dtfrj 



r\' ,a — out 



= W(^)[i-A(^)] 

r),ct—in 

-R V aCME)[l-Me-«)] (9) 
= R^[J V (E) [1 -/„,(£_„)] - 

[!-/,(£)] (10) 



where e± a = E ± hw a and R^ a (^6cJ denotes the 
emission (absorption) scattering rate for the a-phonon. 
Equation[9](Eq.[T0|) is the collision integral for a-phonons 
where the emission process scatters an electron from the 
carrier branch 7/ (77) to the carrier branch 77 (7/). These 
rates have different functional forms depending on the 
phonon branch involved in individual scattering events. 
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FIG. 2: Phonon scattering mechanisms: (a) interbranch AP, 
(b) interbranch OP and (c) intrabranch OP. 



Here we point out that any elastic backscattering, 
e.g. scattering with impurities, would have given similar 
expressions with their specific rates, provided that they 
do not depend on carrier energies. 

The OP emission and absorption scattering rates are 
written as: 



J- l^arm / 



N a (q) + 1 



(17) 



The total contribution to carrier scattering with phonons 
is: 



and 



dtf v 



ph 



E 9 ^ 

a>,r}' 



7]' ,a — in 



£ 

a' ,rj' 



dtfn 



rj' ,a' — out 



(ii) 



For any given phonon wave-vector q and energy 
ftiu a (q), the initial and final carrier wave- vector are de- 
termined by using momentum and energy conservation 
during collision. For instance with the Fermi energy lo- 
cated at the Dirac point (k = in Eq.[Tj), the conservation 
relations are written as: 



E v (k) 



k'±q 

E v ,{k')±huj a {q) 



(12) 
(13) 



where the + (— ) sign denotes phonon emission (absorp- 
tion) process from the fc-state in the branch r] to the k'- 
state in the branch r\ . In our model we consider AP scat- 
tering, that dominates transport in the low-bias regime, 
and OP scattering that induces the electric transport 
nonlinearities in high-field. The different mechanisms are 
depicted in Fig. [2l 

The AP populations are considered to be in thermal 
equilibrium with the lattice as their energies are smaller 
than thermal fluctuations at room temperature. In ad- 
dition, we assume that carrier collisions with APs are 
elastic (i.e. no energy is exchanged with the lattice but 
only between carrier populations) due to the small ratio 
between the sound and Fermi velocity (v s /vf < 1/40). 
Under these assumptions, the AP contribution to the col- 
lision integral vanishes when the initial and final states 
have the same Fermi velocity sign, and emission and ab- 
sorption rates become equal. Therefore, the collision in- 
tegral due to interbranch AP scattering is given by: 



d t f± 



AP 



J2RAp[f T (E)-f ± (E)], (14) 



with Rap = {Z\kpT) / (pVgh 2 vp) by using the deforma- 
tion potential approximation with Za ~ 5 eV— . The first 
two moments of the collision integral for AP scattering 
read—: 



J~ o,± — 



9c 



ttHvf 



£< 



AP 



T)' 



TAP _ 9c^^b ( ( 



2 Tl)J2Ri P . 



(15) 
(16) 



T>V,V _ 
n ab,a ~ 



N a (q) 



(18) 



respectively. In Eqs. fTTl and fT8| r^ p and N a (q) are the 
bare relaxation time and the occupation number for the 
a-phonon, respectively. The experimental estimate for 
the former parameters is about 30 fs&£ and is an order 
of magnitude smaller than those values obtained by first 
principle calculations^. We distinguish between inter- 
branch (rj r]') and intrabranch (77 = 7/) OP scatter- 
ing, as depicted in Fig. O and consider two interbranch 
(ZB and OPi) and one intrabranch (OP2) phonon modes 
interacting with carriers in m-CNTs^. Their energies 
are assumed to be constant and given by Hujzb =0.16eV 
(zone-boundary phonon) and hwop 1 = hcoop 2 =0.20eV 
(interbranch and intrabranch optical phonon), respec- 
tively. 

If the OPs are in thermal equilibrium with the lattice, 
N a (q) = N eq = l/[exp(hcj a (q)/k B T) - 1] as given by the 
Bose-Einstein distribution, where T = T(x) is the local 
lattice temperature. Otherwise, the occupation number 
N a (q) depends on the local temperatures of carriers and 
lattice as well as the current level. At steady state, the 
occupation number of the a-phonon is obtained by using 
the following rate equation: 



d t N a (q) 



ep 



d t N a (q) 



I PP 







(19) 



In Eq.[19]we neglect phonon diffusion because of the small 
OP group velocity. The first term in Eq. [19] accounts for 
the generation/annihilation of OP due to e-p scattering 
and is given by: 



d t N a (q) 



-R^1U(E)[1-Ue^)] (20) 



where E is determined by using Eqs. [12] and [13] The 
second term in Eq. [19] describes the decay of OPs into 
APs and is expressed as: 



d t N a (q) 



-—[N a (q)-N eq (q)} 



(21) 



where OP decay time r pp has been estimated to be 1- 
7ps from Raman studies in m-CNT a 2Q i 27 ^. In addition, 
decay times have been found to vary as the inverse of 

y21,27 
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For intrabranch (77 = rj f ) OP scattering, the phonon 
wave- vectors are q = ±cuop 2 / v f, independently of the 
initial or final state. The moments of the collision integral 
become: 



-rOP 2 
J~ i,± 



-rOP 2 
r 0,± 

= g 2 c (hjop 2 f 
7rhvFTop 2 

NopM 







1 exp 



' k B T ± J 



1 



(22) 



(23) 



For interbranch scattering with ZB and OPi modes, the 
phonon wave-vector q depends on the initial and final 
electron states. Consequently, N a (q) is obtained from 
Eq. [19l and depends on the carrier and lattice temper- 
atures as well as the quasi-Fermi level separation (cur- 
rent). For r pp 7^ and T + 7^ T_, the moments of the 
collision integrals (Eq. [6|) for interbranch OP scattering 
do not have an analytical expression and are obtained 
numerically. In the Appendix we derive the expressions 
for the moments of the collision integral corresponding 
to interbranch scattering with OPs when r pp = and 
T + =T_. 



C. Electron-phonon heat exchange 
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FIG. 3: Flowchart of the heat transfer. The power given 
to the system by the external field is transferred to the 
forward and backward carrier populations. Simultaneously, 
these populations exchange heat with OP populations by 
emission/absorption processes. In addition OPs decay into 
APs, transferring heat to the lattice. The heat transferred 
from the field is removed through the leads or the substrate. 



Combining the 1st moment equation for forward and 
backward carrier populations (Eq. [5]), the heat produc- 
tion/dissipation satisfies: 



Qel + Qlat 



= IF 



where 



qel 



±g c Gj h d x T + +g c G- h d x T- 



(24) 



(25) 



is the amount of heat carried by the electrons and 



Qlat 



hw a {q)d t N a {q) dq (26) 



ep 



is the amount of energy removed from the carriers for 
the OP populations. In addition, since heat transport 
is neglected in Eq. [2lJ qi at denotes the rate (per unit 
length) of energy transferred to the lattice (APs) by OP 
decays. This value determines the local lattice heating, 
and depends on the current and the carrier and lattice 
temperatures. We assume that lattice (AP) heat trans- 
port is diffusive and follows the Fourier's law for which 
the heat equation reads: 



Qap + Qsub — Qlat- 



(27) 



The first term in the left hand side of Eq. [57| accounts for 
the heat carried by the AP population, 



q ap = -Ad x [n(T)d x T] , 



(28) 



The factor A = irdt is the cross sectional area, where d 
and t w 0.34nm are the diameter and the thickness of 
the CNT, respectively. The lattice thermal conductiv- 
ity k(T) is given by k(T) = kqTq/T due to the Umk- 
lapp phonon-phonon scattering 29 . The room tempera- 
ture thermal conductivity has been estimated to be be- 
tween 15 and 60 W/(cm K) for SWCNTs 30 . The second 
term in Eq. [27l determines the heat removal through the 
substrate, which is modeled as a contact resistance 31 : 



Qsub = go(T - T ) . 



(29) 



The parameter go denotes the thermal coupling between 
the CNT and the substrate which has been experi- 
mentally estimated between 0.05 and 0.20 W/(Km) for 
Si02 13,14 . The lattice temperature boundary conditions 
between CNT and leads are set by the presence of a con- 
tact resistance: 



nd x T\ 



T(±L/2) - T 



(30) 



b=±L / 2 nt h 

with Kth ^ 10 7 K/W for m-CNTs^ 2 -. 

D. Heat flow diagram 

The coupled nonlinear integro-differential equations 
for the carrier and lattice temperatures requires a self- 
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consistent solution of Eqs. [5] and [271 The interplay be- 
tween carrier and phonon heat transfer (Eq. is de- 
picted in Fig. El where both forward and backward car- 
rier populations gain energy from the electric field. In the 
high-field regime, OP population builds up (hot phonon 
effects) due to the enhanced scattering by phonon emis- 
sion, but because of their small group velocity, heat re- 
moval by OPs is not significant. Consequently, either 
they transfer their energy back to the carrier popula- 
tion in reabsorption processes, or they decay into APs, 
thereby heating the lattice. Energy can be removed from 
the system through the leads (by both carriers and APs) 
or by the substrate. A key factor that characterizes heat 
transfer is the ratio r^ p /r pp which determines the frac- 
tion of heat transferred from the lattic o 33 i 34 . On the one 
hand if T^ p /r pp ^> 1, most of the energy gained by the 
OP population is transferred to the APs by instantaneous 
phonon decay and carried by the lattice. On the other 
hand, if r^ p /r pp <C 1, the energy of an OP emitted by 
carrier scattering is more likely to return to the carrier 
population by reabsorption rather than decay into APs. 
Hence, most of the energy remains within the carrier and 
OP populations. However, since OPs cannot carry signif- 
icant amount of heat (because of their small group veloc- 
ity) , heat is removed by the electron populations through 
the leads. In this case, the presence of a substrate would 
make no difference - considering that r^ p or r pp are not 
affected by the presence of the substrate - because only 
a small fraction of heat is transferred to the lattice. 



III. RESULTS 
A. Low-field regime 



in m-CNT obtained experimentally^ for the parameters 
^ap,o = 800nm, tzb = 60fs, top 1 = lOOfs Hujzd= 0.16eV, 
huj p 1 = 0.20eV, and T = 300K as shown in Fig.H 
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T [K] 

FIG. 4: Temperature dependence of the low- field resistivity: 
theory (solid lines) and experiment (symbols) from Purewal 
et alM. 

In Fig. [5] we plot the forward carrier (solid line) and 
lattice (dashed line) temperature profiles corresponding 
to the low- field regime (/ = 0.01/iA) for 0.3 and 3/im-long 
m-CNTs. These profiles are determined for both SS and 
FSCNT configuration, for which the thermal coupling to 
the substrate are go ~ 1.5 x 10 _3 W/(cm K) and 0, re- 
spectively. In the case of symmetric boundary conditions 
for carrier temperatures, i.e. T + (— L/2) = T_(L/2), and 



Heating of carrier or lattice population is negligible 
(\T±(x) - T |, \T(x) - T | <C T ) in the low-field regime. 
In this case, the collision integrals for OP scattering have 
an analytic form (see Appendix). Combining Eqs. [4] with 
[T5l and lA.li we obtain the Matthiessen rule for the low- 
field resistivity: 

m =H-k+£) (3i) 

where A^p and A^p are the effective AP and OP mean 
free path given by: 

x% = nx^o (32) 

and 
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foo\ FIG. 5: Forward carrier (solid line) and lattice (dashed line) 

°P J Z-^/ vpT^ sinh [hco a / (2kBT)] 2 j temperature difference to the contact temperature To corre- 

a sponding for / = 0.01 fiA (low- field regime): (a) 0.3 /xm-long 

This analytical expression is in very good agreement with SSCNT, (b) 0.3 /xm-long FSCNT, (c) 3 /xm-long SSCNT, and 

the temperature dependence of the low-field resistivity ( d ) 3 ^ m ~ lon ^ FSCNT. 
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lattice temperature, i.e. T(—L/2) = T(L/2), the carrier 
temperature profiles satisfy T+(x) = T_(— x) due to the 
electron-hole symmetry, thereby reducing the number of 
differential equations. In this regime the electric field 
along the channel remains constant. For short CNTs the 
profiles are similar regardless of the presence/absence of 
substrate with the heating in the carrier population being 
more significant than in the lattice. The forward carrier 
temperature increases along the tube, and the maximum 
is reached at the drain as a consequence of the long mean 
free path (compared to the CNT length)^. In the case of 
long m-CNTs, despite of substantial inelastic scattering 
that heats the AP population, the carrier temperature 
increase is larger than that of the lattice. In particular, 
for FSCNT lattice heating is comparable to that of the 
carriers, while in the SSCNT the lattice temperature is 
flattened by the efficient heat removal through the sub- 
strate, and remains closer to the substrate temperature. 
The maximum lattice temperature is located exactly at 
the m-CNT mid-length, while the forward carrier tem- 
perature peak is shifted towards the right lead. 



B. High- field regime 



t T+ [K] 




0.2 0.4 0.6 0.8 
V ds [V] 



FIG. 6: IV characteristics of 0.3 and 3/im-long CNTs in both 
SS and FS configurations. Symbols correspond to experimen- 
tal data from Javey et al. 9 (□) and Pop et al.— (CO- 

In Fig. [6] we display the electrical characteristics (solid 
lines) corresponding to the SS and FS configurations of 
0.3 and 3/im-long CNTs that are in excellent agreement 
with experiments^^. The parameters used in the cal- 
culations are: ^(300K) = 20 W/(cm K), g = 1.5 x 
10 _4 W/(cm K) and Za ~ 5 eV. The relaxation times 
t zb ~ 85fs, rop 1 = tqp 2 ~ 200fs and r pp = lps are within 
the range estimated by experiments^^ and theory^. 
Nonlinear electrical characteristics are obtained for all 
CNTs. In the high-bias regime, the CNTs in FS configu- 
ration exhibit negative differential resistance (NDR) in- 
duced by the inefficient heat removal that increases scat- 
tering. However, while in 0.3 /im-long FS and SS CNTs 
the current levels approach the 20 /iA limit attributed to 
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FIG. 7: Lattice (left) and carrier (right) temperature profile 
as a function of the bias voltage Vds for the IV characteris- 
tics shown in Fig. [6] From top to bottom: 0. 3/im-long SS- 
CNT, 0.3/im-long FSCNT, 3/im-long SSCNT, and 3/im-long 
FSCNT. 



the onset of OP emission 7 , in the 3 /im-long FSCNT the 
current level is substantially reduced (but still exhibiting 
NDR features). By contrast in the 3 /im-long SSCNT, the 
current increases monotonically in the bias range shown, 
and no saturation is observed due to the fact that the 
electric field is about an order of magnitude smaller than 
in the 0.3 /im-long SSCNT for the same voltages. 

The corresponding lattice and forward carrier temper- 
ature profiles are shown in Fig.0 For short (0. 3/im-long) 
m-CNTs, the temperature profiles in both FS and SS con- 
figurations have a parabolic shape, as shown in Fig. [71(a) 
and 0(b). Unlike the low- field regime (Fig. [5]), the tem- 
perature maxima for forward carriers in the high-field 
regime occur at the center of the m-CNT in both cases, in 
agreement with breakdown experiments in m-CNTs—^. 
We note that for the same channel length, the values of 
the carrier and lattice temperature maxima are lower in 
the SSCNT than in the FSCNT as heat is removed by 
the substrate. In SSCNTs, the reduced temperature low- 
ers scattering rates, and consequently higher current lev- 
els are achieved compared to their FS counterparts. By 
contrast, in 3/im-long SSCNTs the efficient heat removal 
through the substrate flattens the temperature profiles as 
depicted in Fig. [3(c), as observed recently^ 8 . Because of 
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the electron-hole symmetry (T + (x) = T_(— x)) forward 
and backward carrier temperatures are equal to one an- 
other in the channel, but are different in the regions close 
to the contacts. However due to the effective heat trans- 
fer to the substrate, carrier populations are not locally in 
thermal equilibrium with the lattice. For the 3 /im-long 
FSCNT (Fig. [7|d), carrier and lattice temperature pro- 
files have a parabolic shape and markedly higher values 
due to limited heat removal, which enhances carrier scat- 
tering and phonon decays. The current level reduction as 
well as the near thermal equilibrium between carrier and 
phonon populations indicates that transport approaches 
the diffusive regime limit 39 . 



C. Heat flow 

In Fig. [8] we simultaneously display the local heat pro- 
duction given by Joule's law (Eq. [24j) and the fraction of 
the local heat carried by electrons (q e i in Eq. |25|) . APs 
(q ap in Eq. (27j), and substrate (q S ub in Eq. [29]) for the 
IV characteristics shown in Fig. [6j For illustrative pur- 
poses, the profiles have been normalized to the maximum 
power per unit length across the channel for Vd s = IV. 
Since the current remains constant along the channel the 
g-profile is proportional to the electric field profile and 
presents a peak at the middle of the m-CNT. Despite 
significant carrier heating (left column of Fig. [7j), the 
fraction of heat carried by electrons (second column in 
Fig. [8]) is small (< 10%) compared to the total heat pro- 
duced locally for all cases shown. Indeed, most of the 
heat is removed by APs (lattice) and/or the substrate 
(in the SS configuration), as depicted in the third and 
fourth columns of Fig. [HJ For 0.3 /im-long CNTs, the 
heat production profiles in the SS and FS configurations 
look alike despite the different heat removal mechanisms, 
as shown in Fig. [8ja and Fig. [8jb. In contrast for 3 /im- 
long m-CNTs, the field is uniform in the SS configura- 
tion (deviations are less 10% for biases shown), whereas 
it is and highly inhomogeneous with its peak value at the 
CNT midlength in the FS configuration. We notice that, 
in sufficiently long SSCNTs where the substrate flattens 
the temperature profile, the local field value does not 
depend on the boundary conditions. In this case, an ap- 
proximate solution can be obtained by solving for the 
values of carrier and lattice temperatures such that their 
corresponding temperature gradients vanish in the cou- 
pled differential equations (Eqs. 14} 151 and |27|). This self- 
consistent analysis of the thermal flow shows that the 
carrier contribution to heat transport in m-CNTs with 
L > 300nm can be neglected and, therefore, validates 
the usual approximatio n 10 ^ 3 qi a t ~ IF. 

We point out that recent work o 4Q i 41 have suggested 
the existence of surface polariton scattering in the SS- 
configuration. This scattering mechanism would indeed 
modify the heat transport scheme (Fig. [3]) by adding 
a new carrier-substrate channel to heat removal (with- 
out the mediation of APs), which could be accounted for 



by including an additional term on the left hand side of 
Eq. [231 However, the strength of this mechanism, which 
decays as a power law with the CNT-substrate distance, 
and possibly varies with gate bias, still requires confir- 
mation. 



D. Hot Phonons 

By solving Eq.[2l]self-consistently with the carrier and 
lattice temperature profiles for a particular voltage we 
compute the OP occupation number N a (q) and deter- 
mine the equivalent of the temperature of the a-phonon 
mode T a (q) as: 



T a (q) 



k B log 



N a (q) 



(34) 



In Fig. [9l we plot the lattice temperature profile T(x) 
(narrow surface plot), the temperature TzB(q,x) corre- 
sponding to the ZB phonon mode, and the difference be- 
tween the phonon mode temperature and the local tem- 
perature profile, i.e. ATzB(q,x) = TzB(q,x) — T{x) for 
V ds values 0.2, 0.4, and IV in the 3 /im-long CNTs for 
which we have previously computed the electrical charac- 
teristics and corresponding temperature profiles (Figs. [6] 
and[7j). On the one hand we observe that the phonon 
temperature profile in SSCNT is uniform along the in- 
direction (except for the regions close to the leads) be- 
cause of the constant carrier and lattice temperature pro- 
files as displayed in Fig. [91(a). For Vds = IV, the high 
current levels (i.e. large quasi-Fermi level separation) fa- 
vor scattering by OP emission 42 and a non-equilibrium 
OP population builds up (due to non-zero r pp ). On the 
other hand, the OP equivalent temperature profiles in 
the 3 /im-long FSCNT are not uniform reaching peak 
values close to 1000 K for Vds = IV. Nonetheless, in 
the FS configuration the deviations ATzb are smaller 
than in the SS one, as shown (in different color scales) 
at the bottom of Fig. (9](a) and Fig. [S^b), respectively. 
The high lattice temperature achieved in FSCNTs in- 
creases both emission and absorption scattering rates, 
thereby damping (non-equilibrium) hot phonon effects 
as carrier and phonon populations reach local thermal 
equilibrium (T±(x) = T(x) = Top(x)). Therefore, in 
the absence of hot phonon effects, our model shows that 
the cause for the onset of NDR at high bias is, as an- 
ticipated by Conwell 4 ^, the existence of the electric field 
inhomogeneities along the heated CNT 42 . 



E. Electrical power vs length 

Electrical breakdown by oxidation in CNTs has been 
estimated to occur at around 1000K 44,45 . On Fig. [10] 
we display the electrical powers P(1000 K) (solid line) 
and P + (1000 K) (dashed line) as a function of the 
CNT length for various values of the thermal coupling 
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-0.5 
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FIG. 8: 2D heat flow profiles as a function of Vds and distance for 0.3 and 3/im-long CNTs, for the IV characteristics and 
temperature profiles depicted in Figs. [6] and [3 respectively. From top to bottom: 0.3/mi-long SSCNT, 0.3/xm-long FSCNT, 3/xm- 
long SSCNT, and 3/xm-long FSCNT. First column corresponds to the Joule's heat production (qtot)- Second, third and fourth 
column display the fraction of heat carried by the electrons (q e i), acoustic phonons (q a p), and substrate (q su b), respectively. 
The profiles are normalized to the maximum power density across the channel at Vds = IV. 



to the substrate go. The parameters P(1000 K) and 
P + (1000 K) correspond to the electrical power (P = IV) 
dissipated along the CNT when the maximum lattice 
temperature and maximum carrier temperature reach 
1000K, respectively. For the FSCNT (g = 0), both 
electrical powers are determined from the self-consistent 
solution of the BTE and heat equation. For SSCNTs 
(go > 0), we compute the power as P = I(IR q + FL) 
by assuming that the electric field is homogeneous along 
the conductor (vanishing thermal gradients) and solv- 
ing for the values of electric field and currents as de- 
scribed in Section Ull CI In the former case, i.e. go = 
(□), P(1000 K) andP + (1000 K) are approximately equal 
(both curves are indistinguishable when L > 1/im), 
and scale inversely proportional to the length. In con- 
trast the dissipated powers in SSCNTs increase linearly 
with length for sufficiently long m-CNTs, when the volt- 
age drop at the contacts is smaller than that along the 
channel. Similar trends have been reported experimen- 
tally in the electrical breakdown of multi-walled CNTs 
in both FS and SS configuration^. We point out that 
at a fixed length P + (1000 K) and P(1000 K) increase 
with the substrate coupling. Additionally, the difference 
AP = P(1000 K) - P + (1000 K) widens by the enhance- 
ment of go indicating more prominent hot electron effects 
as carrier populations reach the 1000 K temperature at 



lower power values than the lattice. 



IV. CONCLUSIONS 

A model for electro-thermal transport in ID metallic 
systems in the high temperature (incoherent) regime has 
been presented. We use the BTE formalism to describe 
the interdependent electron-phonon dynamics in high- 
field. By accounting for the emergence of carrier and 
phonon populations out of equilibrium we quantify the 
heat production and dissipation by each of these popula- 
tions depending on the experimental configuration. The 
phonon decay rate is the bottleneck that regulates heat 
exchange between carriers and lattice. Our model shows 
remarkable agreement with high-field transport experi- 
ments and offers a qualitative interpretation of break- 
down experiments in multi- walled- CNTs. We attribute 
the emergence of negative differential resistance observed 
in suspended CNTs to the inhomogeneous electric field 
and self- heating in the high field regime. 

One of the authors (MK) acknowledges the support 
of the Department of Physics and the Department of 
Electrical and Computer Engineering at the University 
of Illinois at Urbana-Champaign. 
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(a) SSCNT v = o.2v v ds = o.4v v ds = 1 .ov 




x/L x/L 



(b) FSCNT v ds = o.2V v ds = o.4V v ds = i.ov 




x/L x/L x/L 



FIG. 9: 2D contour plot of the non-equilibrium ZB phonon 
equivalent temperature as a function of the phonon wave vec- 
tor q and distance for a 3/xm-long CNT in: (a) SS and (b) 
FS configuration under various external biases: 0.20V (left), 
0.40V (center), and 1.00V (right). For each bias the lattice 
temperature T(x) profile along the channel, the effective tem- 
perature of the first phonon mode Tzb(x, q) and the tempera- 
ture difference ATzb(x, q) between each mode and the lattice 
temperature are shown in the top, middle and bottom color 
plots, respectively. Note that TzB(x,q) are plotted using the 
same color scale for both lengths but ATzs(x,g) has been 
represented with different ones. 




0.3 1 3 10 

L [urn] 



FIG. 10: Electrical powers, P(1000 K) and P+(1000 K), dis- 
sipated in the CNT when length when the maximum lattice 
(solid lines) and carrier (dashed lines) temperatures reach 
1000K, respectively. The symbols indicate the different values 
of the thermal coupling to the substrate. 



APPENDIX: ANALYTICAL EXPRESSIONS FOR 
THE COLLISION INTEGRALS 

For carriers in different branches in thermal equilib- 
rium (i.e. T + = T_ = T e i), and when hot phonon effects 
are neglected (i.e. r ep ^> r pp or low current level), an- 
alytical expressions for the zeroth and first moments of 
the collision integral (Eq. [6]) and the OP generation by 
phonon emission/absorption (Eq. [20]) can be obtained: 



hlTVFTop 

= _ 9^ak B T el {{N ^ q) + [g{Q+) + _ Na{q) + g( _ n+)]} (A 2) 

ZaiTVFTop 



d t N a 



g*fhj a k B T e 



{(N a (q) + 1) [g(Sl+) + - N a (q) [^(-0_) + G(-Sl+)]} (A.3) 



2h7TVFT op 

r 



where we define Q{x) = exp( ,^_ 1 and Q± = • is valid approximation for the long FSCNTs because the 

The first two expressions are valid in the low-bias regime, minor thermal imbalance between OP and AP (lattice) 

where heating in the carrier population and lattice can be populations, 
neglected. In addition, the omission of hot phonon effects 
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